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ABSTRACT 



We performed spectroscopic observations for a large infrared QSO sample with a 
total of 25 objects. The sample was compiled from the QDOT redshift survey, the 1 Jy 
ULIRGs survey and a sample obtained by a cross-correlation study of the IRAS Point 
Source Catalogue with the ROSAT All Sky Survey Catalogue. Statistical analyses of 
the optical spectra show that the vast majority of infrared QSOs have narrow permitted 
emission hues (with FWHM of H/3 less than 4000 kms~^) and more than 60% of them 
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are luminous narrow line Seyfert 1 galaxies. Two of the infrared QSOs are also classified 
as low ionization BAL QSOs. More than 70% of infrared QSOs are moderately or 
extremely strong Fe II emitters. This is the highest percentage of strong Fe II emitters 
in all subclasses of QSO/Seyfert 1 samples. We found that the Fe II to H/3 line ratio is 
significantly correlated with the [O III] A5007 peak and H/3 blueshift. Soft X-ray weak 
infrared QSOs tend to have large blueshifts in permitted emission lines and significant 
Fc II 48, 49 (5100-5400 A) residuals relative to the Boroson & Green Fe II template. If 
the blueshifts in permitted lines are caused by outflows, then they appear to be common 
in infrared QSOs. As the infrared-selected QSO sample includes both luminous narrow 
line Seyfert 1 galaxies and low ionization BAL QSOs, it could be a useful laboratory to 
investigate the evolutionary connection among these objects. 

Subject headings: galaxies: Seyfert — quasars: emission lines — quasars: general 

1. Introduction 

Two of the most important reasons for investigating the ultraluminous IRAS galaxies (ULIRGs) 
are to find the evolutionary connection between circumnuclear massive starbursts and active galactic 
nuclei (AGNs) and to identify the evolution path from galaxy mergers to elliptical galaxies and QSOs 
(see Sanders & Mirabel 1996 for a review). In recent years, significant progress has been achieved 
with both space and ground-based telescopes (Surace et al. 1998; Surace, Sanders, &; Evans 2000; 
Genzel et al. 1998; Farrah et al. 2001 and reference therein). It is now widely accepted that the 
vast majority (> 95%) ULIRGs are strongly interacting and merging galaxies while some of them 
are post-merger galaxies. The AGN phenomenon probably appears at the final merging stage (e.g., 
Clements et al. 1996; Kim, Veilleux, & Sanders 1998; Zheng et al. 1999; Canalizo & Stockton 
2001a; Cui et al. 2001). Spectral analyses for large samples of ULIRGs reveal that the fraction of 
objects with AGN spectral characteristics is about 25-30% while the fraction of QSOs/Seyfert Is is 
less than 10% (Wu et al. 1998; Lawrence et al. 1999). However, the percentage of QSOs/Seyfert Is 
increases with increasing infrared luminosity, reaching 30-50% for Lir > lO^^'^L© (Veilleux, Kim, 
& Sanders 1999). As the infrared luminosity of ULIRGs is equivalent to the bolometric luminosity 
of optically selected QSOs (Sanders 2001), we find it convenient to refer to QSOs/Seyfert Is selected 
from ULIRGs as IR QSOs throughout this paper. 

Previous spectroscopic studies of small and statistically incomplete IR QSO samples have 
uncovered some unusual properties compared with optically selected QSOs/Seyfert Is. Many 
IR QSOs are extremely strong Fe II emitters (Fe IIA4570/H/3 > 2.0), for example, PEL 1092, 
IRAS 07598-F6508 and Mrk 231 (Li'pari 1994; Lawrence et al. 1997). In fact, almost 100% of 
extremely strong Fe II emitters are luminoTis IR QSOs (Lipari et al. 2002). More than twenty years 
after their discovery, the origin of such extremely strong optical Fe II emissions in QSOs/Seyfert Is 
is still being debated. It has become clear, however, that the strength of the Fe II emission cannot 
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be explained in the framework of photoionization excitation. If strong outflows and shocks are 
present in strong/extremely strong Fe II emitters, non-radiative shock heating and overabundance 
of iron may help to explain the strong Fe II emission (Collin &: Joly 2000). Studies of IR QSOs 
may therefore shed new light on the origin of the Fe II emission. 

Furthermore, the fraction of low ionization broad absorption line QSOs (lo-BAL QSOs) is much 
higher in an IR QSO sample (27%) than that in an optically selected QSO sample (1.4%) (Boroson 
k, Meyers 1992). The lo-BAL QSOs are defined as a subclass of broad absorption line QSOs with 
an obvious Mg II AA2795,2802 doublet and low-ionization line absorption troughs (Weymann et 
al. 1991). Some IR QSOs belong to yet another class of unusual AGNs, the luminous narrow line 
Seyfert Is (hereafter NLSls, Osterbrock & Pogge 1985). NLSls are defined by their optical emission 
line properties. They have narrow hydrogen Balmer lines with typical full width at half maximum 
(FWHM) K, 500-2000 kms'^ The [O III] A5007/H/3 line ratio is less than 3 and most of them have 
strong Fe II emissions. In the X-ray band, NLSls show systematically steeper power-law slopes 
in the continuum than normal Seyfert Is. Some NLSls exhibit rapid X-ray variabilities as well 
(Boiler, Brandt & Fink 1996). Moran, Halpcrn, & Hclfand (1996) pointed out that many NLSls 
are luminous in the infrared band. IR QSOs therefore offer a unique opportunity to investigate 
potential physical connections between IR QSOs, lo-BAL QSOs and luminous NLSls (Brandt & 
Gallagher 2000; Canahzo & Stockton 2001b; Sanders & Mirabel 1996). 

In this paper, we study the spectroscopic properties of IR QSOs based on a large sample of 25 
objects. We compare their properties with those of the Boroson & Green (1992, hereafter BG92) 
sample, which includes 87 optically selected QSOs. The outline of the paper is as follows. In §2, we 
discuss how our IR QSO sample is compiled. The observations and data reduction are described 
in §3. We present the spectra and statistical properties of our sample IR QSOs in §4 and finally, 
in §5, we discuss and summarize our results. Throughout this paper we use a Hubble constant of 
Hq = 50 kms^-*^ Mpc^-*^ and ilg = 1 ^-nd no cosmological constant. As all our objects have rcdshift 
lower than 0.35, the adoption of a different density parameter and cosmological constant has little 
effects on our results. 

2. Sample Selection 

Our IR QSO sample is compiled mainly from three sources 

• The QDOT redshift survey is a survey of the IRAS galaxies sparse-sampled at a rate of one 
in six. It includes 2387 IRAS galaxies complete down to a 60 ^m. flux density limit of 0.6 Jy. 
Lawrence et al. (1999) gave a table of 97 ULIRGs (for Hq = 50 kms^^Mpc-^ and ^^0 = 1) 
with the criterion that the 60 /um luminosity is greater than IO^'^Lq. The optical spectral 
features and classifications are given in the ULIRG table. There are eight objects in the table 
identified as IR QSOs. 
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• Kim k Sanders (1998) selected 118 ULIRGs (infrared luminosities L(8-1000 /xm) > 10^^ L© 

for Ho = 75 kms~^ Mpc~^ and JIq = 0) from the criterion of 60 /im flux density greater than 
1 Jy in sky region oi S > —40° and \b\ > 30°. Veilleux et al. (1999) gave the spectroscopic 
features and classification of 108 out of these 118 ULIRGs. There are 10 IR QSOs among 
these 108 ULIRGs. 

• Moran et al. (1996) presented spectroscopic classifications for a catalogue of IRAS galaxies 
selected from the cross-correlation of the IRAS Point Source Catalogue with the ROSAT All 
Sky Survey by Boiler et al. (1992). This catalogue consists of 241 objects and 80 of them are 
identified as QSOs/Seyfert Is. 11 of these objects are IR QSOs. 

Due to the constraint of the observatory site and instrumental capability, we selected our 
targets by requiring the IR QSOs to be in the northern sky (6 > —30°) and z<0.35 with L(8- 
1000 12m) > 10^2 L0 for Ho = 50 kms"^ Mpc"^ and = 1. In addition, we included F09427+1929 
(Zheng et al. 1999). Taking into account the overlapping sources, the sample consists of 25 
IR QSOs. Table 1 lists the basic parameters. Note that all the infrared luminosities have been 
converted using Hq = 50 kms~^ Mpc~^ and Clo = 1. 

While the sample size is still moderate, it is interesting to put this number in the context 
of the total expected number of IR QSOs in the local universe. The PSCz catalogue provides a 
complete redshift survey of 15411 IRAS galaxies (Saunders et al. 2000). About 900 ULIRGs were 
found, which implies the percentage of ULIRGs in the IRAS galaxy catalogue is about 6%. The 
percentage of IR QSOs among ULIRGs is approximately 10%, as there are 8 IR QSOs among 97 
ULIRGs in the QDOT catalogue and 10 IR QSOs among 108 ULIRGs in the 1 Jy sample. So the 
fraction of IR QSOs in the complete PSCz catalogue is of the order of 0. 6%, i.e., the total IR QSOs 
in the PSCz catalog may be less than 100. Our sample therefore includes roughly one quarter of 
IR QSOs in the local universe. Statistical results based on this quite large IR QSO sample should 
be representative. 



3. Observations and Data Reduction 

Long-slit optical spectroscopic observations were carried out on the 2.16m telescope at the Xin- 
glong station of the National Astronomical Observatories. The observations were mostly performed 
between October 1998 and November 1999 using an OMR spectrograph while some preliminary 
studies were conducted before 1998. For our 1998 and 1999 observations, a Tektronix 1024x1024 
CCD was used giving a wavelength coverage of 4000 A to 9000 A with a grating of 200 A mm . 
The spectral resolution was 9.7 A (2 pixels). The slit width varied from 1.5" to 3.5" to match the 
seeing at the Xinglong station. For the pre- 1998 observations, the instrument setup was slightly 
different. For these observations, a grating of 195 Amm~^ was used and the coverage is from 
3500 A to 8100 A with a resolution of 9.3 A. For F01572-F0009, the spectrum coverage was in the 
range of 3800 ~ 6300 A with a resolution of 4.9 A. The observation log is given in Table 2 listing 
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observation epochs, exposure times, approximate seeings and adopted slit widths. 

Data reduction was performed using IRAF software. CCD reductions included bias subtrac- 
tion, fiat field correction and cosmic-ray removal. Sky light subtractions were accomplished during 

the extraction procedure. Wavelength calibrations were carried out using a He-Ar lamp. The 
resulting wavelength accuracy is better than 1 A. KPNO standard stars were observed for flux 
calibrations. The telluric O2 absorption bands near 6870 A and 7620 A were removed using the 
spectra of the standard stars. 

The IRAF package 'SPLOT' is used to measure isolated emission lines. 'SPECFIT'^, an 
interactive spectral analysis procedure linked to IRAF is used to measure blended lines (e.g.. Ha and 
[N II]AA6548,6583). It can match a wide variety of emission lines, absorption lines, and continuum 
models. We model the emission lines with Gaussian profiles and the local continuum as a power-law. 

Since the intrinsic extinction of ULIRGs is significant, the Galactic extinction is ignored in our 
analysis. The extinction correction is calculated approximately according to Veilleux & Osterbrock 
(1987) and the intrinsic broad line Ha/H/3 ratio is taken to be 3.1 (Baker 1997). The H/3/H7 ratio 
is used to estimate the extinction when Ha is out of the spectral coverage (Osterbrock 1989). The 
measured fluxes of Balmcr lines arc the sum of the narrow and broad components. It is noted 
that the intrinsic broad Balmer line ratios may be larger than the adopted values (e.g. MacAlpine 
1985) and the extinction corrections may have large uncertainties. Fortunately, the emission line 
ratios adopted in our analyses, e.g. Fe IIA4570/H/3 and [O III] A5007/H/?, are almost independent 
of extinction, due to the adjacency of the involved lines. For the Fe II residual measurement, 
the uncertainty from the extinction correction will be discussed in section 4.3. In general, the 
uncertainty for the flux measurement introduced by the extinction correction is less than 20%. 

As the Fe II emission is moderately or extremely strong for most of our sample galaxies, 
Fe II multiplets seriously blend with the H/3 and [O III] AA4959,5007 lines and contaminate the 
continuum. We carefully remove the Fe II multiplets following BG92. Their method uses an Fe II 
template derived empirically from high quality data of I ZW 1, a typical NLSl galaxy. We then 
measure the flux for all emission lines based on the Fe Il-subtr acted continuum. On the other hand, 
the emission lines for a quite large fraction of our sample IR QSOs show a remarkable asymmetric 
profile. In such cases, double-Gaussian profiles are needed to fit them. The asymmetry and blueshift 
can then be measured from the fitting. In the following subsections we describe in more detail the 
Fe II multiplets removal method and the measurement of the line asymmetry and blueshift. 



^SPECFIT is developed and kindly provided by Gerard A. Kriss. 
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3.1. Fe II Multiplets Removal 

In order to estimate the Fe II strength and measure the hne fluxes rehably, we adopt the BG92 
method which rehes on an Fe II template. The template and observed spectra are both transformed 
into the rest frame. The template is broadened by convolving with a Gaussian of various line widths 
and scaled by multiplying by a factor indicating the line strength. The best match is then searched 
for in the two-dimensional parameter space of the line width and line strength. A good Fe II 
subtraction is found when the parts of the continuum between the H7 and H/3 and between 5100- 
5400 A (which covers the Fe II multiplets 48, 49) are flat. The best-fit Fe II template emission 
lines are then subtracted from the observed spectrum. The Fe II flux is then determined from the 
best fitting Fe II template between the rest wavelengths 4434 A and 4685 A. The procedure of 
Fe II subtraction is illustrated in Fig. 1. In each panel of Fig. 1, the top curve is the dereddcncd 
spectrum, the bottom curve is the Fe II template while the middle curve shows the Fe Il-subtracted 
spectrum. Note that the Fe Il-subtracted spectra have been shifted downwards for clarity. 

In the BG92 method, it is assumed that the relative strengths of the Fe II lines (within each 
multiplet and among multiplets) are the same for different objects. For half of our spectra, the 
Fe II emission can be subtracted very well by the Fe II template. However, for F00275— 2859, 
IRAS 07598+6508, F09427+1929, Mrk 231 and F20036-1547 , the Fe II multiplets 48,49 (5100- 
5400 A) are stronger than Fe II multiplets 37, 38 relative to the BG92 template. In contrast, for 
F02065+4705, F10026H-4347, Z11598-0112, F20520-2329 and F22454-1744, the Fe II multiplets 
37, 38 (4500-4680 A) are stronger than Fe II multiplets 48, 49 relative to EG template. For these 
objects, significant Fe II residuals can be seen in their Fe Il-subtracted spectra. The extreme 
case is Z11598-0I12, for which nearly half of the Fe IIA4570 (i.e. Fe II multiplets 37,38) is 
left after the multiplets 48,49 are subtracted. Fig. 2 shows two examples of such remarkable 
Fe II emission residuals in the Fe II multiplets 48,49 and 37,38 in the top and bottom panels, 
respectively. Compared with the optical QSO sample of BG92, the IR QSO sample contains more 
objects showing large deviations from the Fe II template. 

The spectra For F13218-h0552 and F23411-F0228 arc too noisy to detect the Fe II lines rehably, 
so we take the Fe II strength for F13218-h0552 from RemiUard et al. (1993) and ignore F23411-I-0228 
in our statistical analysis concerning the Fe II strength. 

3.2. Emission Line Fitting and Measurement 

As a first step, we use a single Gaussian profile to fit each emission line for all target galaxies. 
It works well for most emission lines. However, a single Gaussian profile cannot fit some lines, e.g., 
those with asymmetric profiles. In such cases, two Gaussian components are used to fit the emission 
lines, namely, one narrow component and one centroid-shifted broad component. The blueshift (or 
redshift) is defined as the shift of the broad component relative to the narrow component in units 
of kms~^. Fig. 3 illustrates the two Gaussian component fitting for the permitted line H/3 (top 
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panel) and forbidden line [O III] A5007 (bottom panel). In Fig. 3, the observed and fitted profiles 
are shown by the solid and dashed lines, respectively, while the dot-dashed line is for the Gaussian 
components and the dotted line represents the fitting residual. As one can see, for F01572-I-0009, 
the blueshift of the permitted emission line H/3 of F00275— 2859 is 750 kms~^ while the blueshift 
of [O III] A5007 is 510 kms'^ 

Moreover, there is no correlation between the FWHM of the narrow Gaussian component of 
the permitted line H/3 and that of the forbidden line [O III] A5007. It implies that the narrow 
Gaussian component of the permitted emission lines is different from the narrow lines (such as 
[O III] A5007) from the narrow line region. Efforts were also made to separate the narrow Gaussian 
component contributed by the narrow line region. However, the contribution from the narrow line 
region is usually less than 3% so we ignore this component in our discussions. 

We also measure the asymmetry parameter defined by de Robertis (1985): 

_ Ae(3/4) - Ae(l/4) 
AA(l/2) ' 

where Ac(l/4) and Ac(3/4) are the wavelength centers at 1/4 and 3/4 of the maximum, respectively, 
and AA(l/2) is the FWHM. The asymmetry parameter is positive (negative) if there is excess light 
in the blue (red) wing. 

Based on the Fe Il-subtracted spectra and Gaussian fitting of emission lines, we measured the 
fiux, FWHM, and equivalent width (EW) for each strong emission line for all our objects. The 
fiuxes of broad lines or asymmetric lines refer to the sum of double components. For most targets, 
the uncertainty of flux measurement of emission lines is about 10%, but for low S/N cases it could 
be up to 20%. The uncertainty of the blueshift measurement is within 150 kms~^. However, for 
three sources (F02054+2835, F13218 +0552 and F23411+0228), the fluxes of emission lines have 
large uncertainties due to the poor S/N and hence they are less reliable. 

The dereddened and Fe Il-subtracted spectrum for each sample IR QSO is shown in Fig. 1. 
Table 3 lists the FWHM, blueshifts and asymmetry parameters for H^ and [O III] A5007. The 
intrinsic FWHM values arc obtained from the subtraction, in quadrature, of the observed FWHM 
and that of the instrumental profile, measured from the comparison lamp lines. Note that only the 
significant [O III] A5007 blueshifts (the [O III] A5007 blueshift > 500 kms"^) are listed. In Table 4, 
color excess E(B— V), the equivalent widths of the emission lines and various line ratios are listed. 



4. Statistical Results 

As discussed in the introduction, our sample is a unique one to investigate the physical con- 
nection among IR QSO, luminous NLSls and lo-BAL QSOs. We therefore performed statistical 
studies similar to those in BG92 for their optical sample, which includes 87 QSOs from the Bright 
Quasar Survey Catalogue with redshift less than 0.5 (Schmidt & Green 1983). The results for these 
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two samples will be compared in order to understand any possible evolutionary connection between 
IR QSOs and classical QSOs. 

4.1. The Percentage of NLSls and Strong Fe II Emitters 

It is obvious from Table 3 that for all our IR QSOs except F16136+6550 and F18216+6419, 
the H/3 FWHM is less than 4000 kms~^. This differentiates IR QSOs from classical QSOs as the 
main characteristic of classical QSOs is the presence of broad permitted emission lines with typical 
FWHM between 4000 and 10000 kms"^ (Rodriguez- Ardila et al. 2000). To make it clear, Fig. 4 
shows the distribution of the H/3 FWHM for our sample (top) and the BG92 sample (bottom). 
A Kolmogorov-Smirnov (K-S) test indicates a probability of 9.5 x 10~^ for the two distributions 
being the same. From Table 3 and Fig. 4, the percentage of IR QSOs with H/? FWHM less than 
2000 kms-i in our sample is 60% (15/25). In comparison, only 23% of BG92 QSOs have FWHM 
H/3 less than 2000 kms^^ and one third of BG92 QSOs are classical QSOs with H/3 FWHM larger 
than 4000 kms^^ 

As part of our IR QSOs are selected from the IRAS-ROSAT cross-correlation catalogue, our 
sample may be biased to include more NLSls (Stephens 1989). To check this, we performed 
statistics for the sub-sample of 15 objects selected from two purely IR-selected samples (the QDOT 
redshift survey and 1 Jy ULIRG sample). Out of these 15 objects, 8 are identified as NLSls (53%). 
Therefore, there is no clear difference in the fraction of NLSls between the whole sample and the 
purely IR-selected sub-sample. 

We also carefully investigated whether the H/3 blueshift could influence the H/3 FWHM and 
concluded that this possibility is unlikely. Therefore the percentage of NLSls seems genuinely high 
in our IR QSO sample. 

Fig. 5 shows the histograms of the Fe IIA4570/H/3 ratio for the IR QSO and BG92 samples. 
We can see from Fig. 5 that the Fe IIA4570/H/3 ratio distributions for the two samples are quite 
different. A K-S test reveals that the probability for these two distributions being the same is 
3.0 X 10-^. Table 4 shows that 25% (6 in 24) and 46% (11 in 24) of IR QSOs are respectively 
extremely strong (Fe IIA4570/H/3 > 2.0) and moderately strong (1.0 < Fe IIA4570/H/3 < 2.0) 
Fe II emitters (Joly 1991; Veron-Cetty, Veron-Cetty, & Concalves 2001). In contrast, for the 
BG92 sample, only 15 out of 87 (17%) QSOs are moderately strong Fe II emitters and there are 
no extremely strong Fe II emitters present. The high percentage of strong Fe II emitters in our 
sample can also be seen by comparing with that of the overall AGN population. Veron-Cetty et 
al. (2001) found that moderately strong Fe II emission occurs in only about 5% of AGNs and 
only a few AGNs are extremely strong Fe II emitters. Furthermore, the two weakest Fe II emitters 
(with Fe IIA4570/H/3 < 0.5) in our sample (F16136+6550 and F18216-F6419) have broad permitted 
emission lines with FWHM larger than 4000 kms~^, and hence are classical QSOs. 
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4.2. The Correlations 

In this subsection, we study the correlations between various emission hnes and continua for our 
sample IR QSOs. For this purpose, we performed Spearman Rank-order (S-R) correlation analyses 
among various quantities and investigate the implications of these correlations. As mentioned 

above, most of our sample IR QSOs arc moderately /or extremely strong Fe II emitters, and the 
origin of the Fe II emission is still not understood. Hence our analysis focuses on the correlations of 
the Fe II strength with other parameters. Throughout this paper, the correlations are characterized 
by the probability V that the null hypothesis of no correlation is true. 

Fig. 6 shows the Fe IIA4570/H/3 line ratio versus the [O III] A5007 peak, defined as the peak 
height of the [O III] A5007 line relative to that of H/3 as in PG92. These two parameters are anti- 
correlated with a correlation coefficient of 0.71 at a very high significance level {V = 6.1 x 10~^). 
This result is consistent with BG92 although the correlation between these two parameters is 
stronger than the one reported by BG92 for optically selected QSOs. The Fe IIA4570/H/3 ratio is 
also well correlated with the [O III] A5007/H/? ratio with P = 9.8 x 10"^. Note that 6 IR QSOs in 
our sample have very weak [O III] A5007 emission (with EW of [O III] A5007 less than 5 A). Such 
a weak [O III] A5007 emission is an important characteristic of lo-BAL QSOs (Boroson &; Meyers 
1992); we return to this point in §5.2. 

Fig. 7 shows the striking correlation between the Fe II A4570/H/3 line ratio and the H/3 blueshift. 
The S-R correlation coefficient is 0.75 which is highly significant with 'P = 2.7 x 10~^. In particular, 
it can be seen from Table 3 and Fig. 7 that about two thirds of our targets show H/? blueshifts with 
values as large as 2000 kms~^. 

As there are striking correlations between both the Fe II A4570/H/3 and [O III] A5007 peak and 
Fe II A4570/H/3 and H/3 blueshift, we also performed a correlation analysis between the [O III] A5007 
peak and H/3 blueshift. We found that these two parameters are well correlated with an S-R 
correlation coefficient of 0.72 corresponding to "P = 1.4 x 10~^. These tight correlations between 
the H/3 blueshift, Fe II A4570/H/3 and the [O III] A5007 peak must reflect some physical connection 
between these lines. 

Fig. 8 plots the Fe II A4570/H/3 line ratio versus the H/3 FWHM for 24 of our objects (F23411+0228 
is excluded due to its low S/N). If the two classical QSOs with H/3 FWHM larger than 4000 km s~^are 
excluded, then the S-R correlation coefficient between these two parameters is 0.47 with V = 
2.8 X 10-2. If we further exclude the 3 objects with H/3 FWHM larger than 3000 kms-^(but 
smaller than 4000 kms"^), the correlation becomes somewhat stronger, with an S-R correlation 
coefficient of 0.58 corresponding to V = 9.9 x 10~^. Wc caution, however, that if we include the 2 
objects with H/3 FWHM larger than 4000 kms"\ the Fe II A4570/H/3 line ratio and the H/3 FWHM 
are no longer well correlated. 

Fig. 9 shows the Fe IIA4570/H/3 line ratio versus the H/3 asymmetry parameter defined by dc 
Robertis (1985, see also eq. 1). The S-R correlation coefficient is 0.46 with V = 2.8 x 10^^. Our 
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result is broadly consistent with BG92, although the H/3 asymmetry parameters for our sample IR 
QSOs seems to be much larger than those for the BG92 sample. Quantitatively, just 13% of BG92 
QSOs have the asymmetry parameter larger than 0.1, while in our sample 48% (26%) of IR QSOs 
have the H/3 asymmetry value larger than 0.1 (0.2). 

As described in §3.2, the H/3 blucshift value is determined by the blueshift of the broad Gaussian 
component relative to the narrow Gaussian component in permitted emission lines. The blueshifted 
broad Gaussian component may be connected with the outflow of clouds in the broad line region 
(Leighly 2001). Outflows with a velocity of several hundred or even a few thousand kms^^ could 
produce shocks that can excite the emission lines. Such shocks may be an important ingredient for 
understanding the strong correlation seen between the Fe II A4570/H/3 ratio and the H/3 blueshift. 

4.3. The Soft X-ray Properties 

We collected all the X-ray information from Moran et al. (1996) and Xia ct al. (2001) based on 
ROSAT archive data. Table 1 lists the ratio of the soft X-ray luminosity to the infrared luminosity, 
Lx/LpiR. This is an important quantity as nearly all the luminosity of an ULIRG is emitted in 
the far-infraxed band (Surace et al. 2000). We can see from Table 1 that 17 of our IR QSOs 
were detected by either ROSAT All Sky Survey (13/17) or pointing observations (4/17). The ratio 
of the soft X-ray luminosity to the far-infrared luminosity spans about four orders of magnitude. 
Not surprisingly the four objects (IRAS 07598+6508, Mrk 231, F00275-2859 and F21219-1757) 
detected only by ROSAT pointings have the lowest Lx/Lfir< 0.01 values compared with the other 
14 objects detected by the ROSAT All Sky Survey. For the non-detections, we calculated their 
upper limits of the soft X-ray luminosity as follows. We assume a detection limit is 6 source 
photons (as in the ROSAT All-Sky Survey Faint Source Catalogue, Voges et al. 2000) and an 
average exposure time of 400 seconds for each source. The spectrum of each source is taken to 
be a power-law with a photon index of 2.3 (appropriate for AGNs) and Galactic neutral hydrogen 
column density is adopted. The estimated Lx/LpiR upper limits are also given in Table 1. As 
might be expected, they all satisfy Lx/Lfir< 0.01. 

Fig. 10 shows the H/3 blueshift versus the ratio of soft X-ray luminosity to far-infrared luminos- 
ity for 14 objects detected by the ROSAT All Sky Survey and ROSAT pointings. In this statistic 
and also in the following statistics concerning Lx/LpiR, we excluded F16136-I-6550, F18216-I-6419 
because they have broad emission lines and are QSOs; F2341 1-1-0228 is also excluded due to its low 
S/N spectrum. It is clear from Fig. 10 that the objects with larger H/3 blueshifts tend to be soft 
X-ray weak and vice versa (the S-R correlation coefficient -0.61 corresponding to V = 2.0 x 10~^). 
In Fig. 10, we also labelled the two potential low-BAL QSOs (F09427-M929 and F20036-1547, see 
section 5.2) using their soft X-ray luminosity upper limits. These two objects strengthen the trend 
described above. 

During our data reduction, we noticed that the Fe II emission cannot be subtracted very well 
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by the Fe II template of BG92 for one third of our targets, i.e., there are large residuals of Fe II 
multiplets 48, 49 (5100-5400 A) or Fc II multiplets 37, 38 (4500 - 4680 A) in the Fe Il-subtracted 
spectra. We define the residual of Fc II multiplets 37, 38 (or Fe II 48, 49) as the excess of Fc II 37, 38 
(or Fe II 48, 49) relative to its strength in the best-fit Fe II template. A positive residual means that 
there is an excess in the Fe II 37, 38 while a negative value signals an Fe II 48, 49 excess. The Fe II 
residual is somewhat uncertain due to the extinction correction. Since the blue light suffers more 
extinction than the red light, an over-correction in the extinction would lead to an over-estimate 
of the Fe II multiplets 37, 38 and hence increase the value of the Fe II residual in this multiplets. 
This overall trend also establishes that the excess of Fe II multiplets 48, 49 relative to multiplets 
37,38 cannot be an artifact. For example, for F00275— 2859, we have adopted E(B— V) = 0.41, 
which results in a negative Fe II residual (i.e., there is an excess in the Fe II multiplets 48,49). If 
no extinction is adopted, then the value of the Fe II residual in the multiplets 48, 49 would further 
increase by 30%. We note that while Boroson & Green (1992) have described such residuals, they 
are not as remarkable as in our IR QSO sample. 

We list the measurements of significant Fe II residuals in Table 4. We have examined the 
correlation of the residuals with other parameters. A correlation was found between the residuals 
of Fe II multiplets and the soft X-ray luminosity. This is interesting because it ties in with the 
long-standing puzzle why some NLSls with extremely strong Fe II emission are X-ray luminous, 
while others with similar optical spectra are X-ray weak. Fig. 11 shows the residuals of Fe II 
multiplets vs. the Lx/LpiR ratio for 14 object with concrete X-ray detections. We can see from 
Fig. 11 that as the residual of Fe II multiplets increases from minus to positive, the soft X-ray 
luminosities of IR QSOs also increase (the S-R correlation coefficient is 0.44 with the correlation 
significant of "P = 0.117). The most interesting result from Fig. 11 is that objects with larger Fc II 
multiplets 48, 49 residuals are all X-ray weak. To check the validity of this result, we measured or 
estimated the Fe II multiplets residual for another two low-redshift lo-BAL QSOs — PG 1700-1-518 
and IRAS 14026-1-4341. The Fe II multiplets 48,49 residuals are 0.19 and 0.1, respectively. These 
two objects are shown in Fig. 11 as well using their upper limits of the soft X-ray luminosity 
(the soft X-ray information of PG 1700-1-518 is from Wang, Brinkmann, & Bergeronet 1996). In 
the same figure, we also indicated the locations of two potential lo-BAL QSOs, F09427-I-1929 and 
IRAS 20036—1547, as they have significant Fe II multiplets 48,49 residuals and extremely weak 
[O III] A5007. It is interesting that all four low-redshift lo-BAL QSO and the two potential low- 
BAL IR QSOs are located in the bottom left of Fig. 11. Furthermore, we also measured the Fe II 
multiplets residual for soft X-ray weak QSOs (Brandt, Laor & Wills 2000) with available spectral 
data. None of them have positive Fc II multiplets residuals. It appears that the Fe II multiplets 
residual may be a good criterion to select X-ray weak QSOs or lo-BAL QSOs. 
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4.4. Infrared Properties 

The infrared color-color diagram has been used as an important tool to discriminate starbursts 
and AGN activities in the nuclear/circumnuclear regions of galaxies (de Grijp et al. 1985; Lfpari 
1994). 20 IR QSOs in our sample are securely detected in three far-infrared bands (25, 60 and 
100/Ltm). Fig. 12 shows the location of these objects in the infrared color-color diagram, a(60, 25) 
versus a(100,60). Here a(Ai, A2) = — log(F(A2)/F(Ai))/log(A2/Ai), where the wavelength is in 
units of iim. In the same diagram, the power-law and blackbody lines are also indicated. 

It is clear from Fig. 12 that almost all IR QSOs except F13218+0552 are located between the 
blackbody and power-law lines. Moreover, there are a group of IR QSOs clearly located close to the 
blackbody line. This group includes F01572+0009, IR06269-0543, F11119+3257, F13218+0552, 
F15462-0450, F23411-F0228 and Mrk 231. We find that they have either significant [O III] A5007 
blueshifts (see Table 3), or significant [O II] A3727 blueshifts (Mrk 231, see Lfpari et al. 2002). 

The objects close to the power-law line are F18216-F6419, F16136-F6550, F12265-F0219 (3C 273), 
F10026-I-4347, F22454-1744 and three additional objects in the bottom left of the figure. The first 
2 objects are the only classical QSOs in our sample while the remaining objects are all moder- 
ate/extremely strong Fe II emitters with bright soft X-ray emission (see Table 4). 

5. Summary and Discussion 

We studied an IR QSO sample with a total of 25 objects. The sample is compiled from the 
QDOT redshift survey, the 1 Jy ULIRGs survey and from a cross-correlation study of the IRAS 
Point Source Catalogue with the ROSAT All Sky Survey Catalogue. Using the observed optical 
spectra and archive data in the infrared and soft X-ray, we investigated the correlations of the 
Fe IIA4570/H/3 ratio with the [O III] A5007 peak, H/3 blueshift and the H/3 FWHM. All these 
parameters are correlated. We found that soft X-ray weak QSOs, especially lo-BAL QSOs, tend 
to have significant Fe II multiplcts 48.49 (5100-5400 A) residuals (cf. Fig. 11). The correlation 
between the Fe II A4570/H/3 residuals and the Lx/Lfir ratio, shown in Fig. 11, although somewhat 
weak, may be a useful clue for understanding why some strong/extremely strong Fe II emitters are 
X-ray luminous while others are X-ray quiet. 

5.1. The Outflows 

One of the striking features of the emission lines for IR QSOs is the blueshift of permitted 
emission lines or the forbidden [O III] A5007 line. As we have removed the Fe II multiplets carefully 

surrounding the H/3 and [O III] AA4959,5007 , there should not be much Fe II multiplet contamina- 
tion to the H/3 and [O III] AA4959,5007 emission lines. Moreover, the blueshifts measured from H/3 
and Ha are similar, as are the blueshifts measured from the [O III] A5007 and [O III] A4959 lines (see 
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Fig. 3). The asymmetries of the H/? and [O III] A5007 hnes cannot be attributed to contaminations 
from other emission hnes. Outflows give a plausible explanation for the emission line blueshift (for 
alternative explanations, see Brandt et al. 2000). Such emission-line outflows have been found 
and discussed extensively for some NLSls (Leighly 2000; Christopoulou et al. 1997), radio galaxies 
(Tadhunter et al. 2001) and for IR QSOs (Lfpari 2002). 

From the strong correlation between the Fe II FWHM and the H/? FWHM, Boroson &: Green 
(1992) suggested that the Fe II line and permitted emission lines share a common emission region. If 
the H/3 blueshifts are due to cloud outflows in the broad line region, shocks are likely to be produced 
in outflows. Such shocks may be responsible for the Fe II emission, as the pure photoionization 
model fails to explain the strong Fe II emission (Collin & Joly 2000) . The tight correlation between 
the H/3 blueshift and the Fe IIA4570/H/3 ratio can then be understood as they are physically 
connected through shocks associated with outflows. Analogously, the [O III] A5007 blueshifts could 
probe the outflows in narrow line emission region. It is interesting to see that all objects with large 
[O III] A5007 blueshifts are located close to the blackbody line in the infrared color-color diagram. 
Hence the locations of IR QSOs in the infrared color-color diagram may be related with outflows 
in the [O III] A5007 emission line region. We caution, however, that there are other explanations 
concerning the strong Fc II emission, especially in NLSl galaxies (see Sulentic, Marziani, & Dultzin- 
Hacyan 2000 for a review). 

If the emission line blueshifts can indeed be attributed to outflows, then our sample indicates 
such outflows are common for IR QSOs. Our statistics do not however clarify which mechanism 
(central AGN radiative pressure, starbursts or both) drives these outflows. High resolution obser- 
vations and investigations in the UV, soft X-ray and optical bands are needed to further explore 
this unique sample. 

5.2. The Connection of IR QSOs with Luminous NLSls and Low- ionization BAL 

QSOs 

60% of our IR QSOs satisfy the strict criteria of NLSls. For the remaining objects, most 
are moderate strong/extremely strong Fe II emitters. However, 71% (5/7) of the extremely strong 
Fe II emitters in IR QSO sample fail to meet the strict criterion of NLSls. Similarly, a large 
fraction (44%) of extremely strong Fe II emitters in the sample of Veron-Cetty et al. (2001) also 
have H/? FWHM larger than 2000 kms""^ As Veron-Cetty et al. (2001) pointed out, there should 
be a continuous distribution of optical line widths for QSOs/Scyfcrt Is and hence the separation 
between the broad line Scyfert Is and NLSls at 2000 kms^^ could be arbitrary (see also Sulentic 
et al. 2000). Our results seem to support the above statement and hint that the strength of 
optical Fe II emission line, rather than the H/3 FWHM, may be the most important characteristic 
of non-classical QSOs. 

Much attention has been paid to lo-BAL QSOs recently. Canalizo & Stockton (2001b) stud- 
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ied four lo-BAL QSOs currently known at z < 0.4 and found that all four are ULIRGs which 

reside in dusty starbursts or post-starbursts. Two of these, Mrk 231 (i.e., F12540+5708) and 
IRAS 07598+6508, are in our IR QSO sample. The main characteristics of lo-BAL QSOs are 
strong Fe II emissions, very weak [O III] A5007 emissions and low X-ray luminosity. They are inter- 
mediately located between the power-law and blackbody lines in the infrared color-color diagram 
(cf. Fig. 12). A recent Chandra survey for BAL QSOs also revealed that all these four lo-BAL 
QSOs have high column densities and are unusually faint in both soft and hard X-ray bands (Gal- 
lagher et al. 2002). These properties are shared by some other IR QSOs in our sample in addition 
to Mrk 231 and IRAS 07598+6508. For example, F00275-2859, F09427+1929 and F20036-1517 
have all the properties of lo-BAL QSOs as described above. We identify them as potential or more 
evolved lo-BAL QSOs. If UV observations establish that they are real lo-BAL QSOs, the percent- 
age of lo-BAL QSOs in our sample will be 20%, similar to the fraction found by Boroson & Meyers 
(1992) for a much smaller IR QSO sample. 

Brandt & Gallagher (2000) argued that the potential physical connection between luminous 
NLSls and lo-BAL QSOs is their high accretion rate relative to the Eddington accretion rate. IR 
QSOs tend to appear in the final merging phase with their central AGN activity recently triggered 
or rejuvenated by the merging activity (Sanders & Mirabel 1996; Zheng et al. 1999; Canalizo &: 
Stockton 2001a and reference therein). Numerical simulations show that a large amount of gas 
flows toward the center during mergers (e.g., Barnes &i Hernquist 1991), so it seems plausible 
that IR QSOs also have high accretion rates. Such systems may have a greater ability to drive 
radiative outflows. The correlations we presented in section 4.2 suggest that the outflow velocity 
and the physical condition in the outflow region (inferred from the H/3 and [O III] A5007 emission 
line blueshifts) influence the Fe II , infrared and soft X-ray emission properties. In short, IR QSOs, 
luminous NLSls and lo-BAL QSOs may have some common physical conditions. The difference 
between them may be from different viewing angles or from different evolution phases. A careful 
study of IR QSOs may be important for understanding the evolution and formation of classical 
QSOs. 

We are grateful to Drs. Th. Boiler, T.Q. Wang, J.Y. Wei, L.C. Deng and Mr. Z.H. Shang for 
helpful discussions, and to the BATC members for advice on data deduction. Particular thanks 
are due to Drs T. Boroson & R. Green for kindly providing their Fe II template and their dataset 
of QSOs in PG92. We also thank the anonymous referee and Dr. Richard James for constructive 
comments that improved the paper. This project was supported by NSF of China and NKBRSF 
G19990754. 
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Fig. 1. — Dereddened spectra (top curves), Fe Il-subtracted spectra (middle curves) and Fe II 
spectra (bottom curves) of 25 IR QSOs, shown in increasing order of right ascension. The Fe II- 
subtracted spectra have been shifted downwards for clarity. The vertical axis shows the observed 
flux in units of 10~^^ ergcm"^ s""*^ A~^. The horizontal axis is wavelength in A in the rest frame. 
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1. — Continued. 
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Fig. 1. — Continued. 
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Fig. 2. — Illustration of significant Fe II multiplets residuals. Top panel: F12540+5708 shows the 
Fe II multiplets 48, 49 residual compared with the Fe II Boroson &: Green (1992) template. Bottom 
panel: Z11598— 0112 shows Fe II multiplets 37,38 residual. In each panel, the top curve is the 
dereddened observed spectrum between 4300-5100 A. The middle one show the Fe Il-subtracted 
spectrum (with an arbitrary offset) with the fitted low-order polynomial continuum shown as a thick 
solid line. The excess of Fe II emission lines is illustrated in the continuum-subtracted bottom curve. 
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Fig. 3.— Illustration of the blueshifted H/3 and [O III] AA4959,5007 lines. Top panel: For 
F00275— 2859, H/3 can be well fitted by a sum of two Gaussian components with the broad com- 
ponent blueshifted by 750 kms~^ relative to the narrow one. Bottom panel: For F01572+0009, 
[O III] AA4959,5007 shows a narrow component in addition to a blueshifted broad component. The 
solid line is the observed profile. The dashed line shows the fitted profile. Each fitted component 
is shown by a dot-dashed line and the dotted line illustrates the residual of fitting. 
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Fig. 4. — Histograms of the H/3 FWHM for IR QSOs (top panel) and for the Boroson & Green 
sample (bottom panel), which includes 87 optically selected bright QSOs. 
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Fig. 5. — Histograms of the Fe IIA4570/H/3 ratio for IR QSOs (top panel) and the BG92 sample 
(bottom panel). 
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Fig. 6.— Fe II A4570/H^ vs. the ratio of the peak height of [O III] A5007 to that of H/?. Upper 
hmits for the non-detections are indicated by arrows. For the detected objects, the Spearman 
Rank-order correlation coefficient is —0.78 with V = 6.1 x 10~^, i.e., the probability for the null 
hypothesis of no correlation between these two parameters being true is 6.1 x 10~®. 



-26- 



o 




-2000 

Blueshift (km s-i 



2000 



Fig. 7. — Fe IIA4570/H/3 vs. the H/3 blueshift, which is defined as the blueshift of the H/3 broad 
component relative to the H/3 narrow component in unit of kms~^. A negative value indicates a 
redshift. For 23 objects, the correlation coefficient is 0.76 with V = 2.7 x 10~^. 
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Fig. 



■ Fe II A4570/H/? vs. H/3 FWHM. The 22 Objects with H/3 FWHM less than 4000 kms^^are 



shown as filled squares while the two objects with H/3 FWHM greater than 4000 kms^^are shown 
with open squares. The Spearman Rank-order correlation coefficient for the objects with FWHM 
smaller than 3000 km s'^ (4000 kms"^) is 0.47 (0.58) corresponding to V = 2.8 x lO'^ (V = 
9.9 X 10^"^). No correlation is apparent if we include all objects. Note that F23411+0228 is 
excluded due to its low S/N spectrum. 
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Fig. 9. — Fe IIA4570/H/3 vs. the H/? asymmetry index defined by de Robertis (1985, cf. eq. 1). 
The correlation coefficient is 0.46 corresponding to V = 2.8 x 10"^. 
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Fig. 10. — H/? blueshift vs. soft X-ray luminosity for 14 objects (filled squares) detected by the 
ROSAT All Sky Survey or Pointings. The correlation coefficient for the detected objects is —0.61 
with V = 2.0 X 10~^. Additionally, two potential lo-BAL QSOs are marked by arrows with their 
upper limits of X-ray luminosity (see text). Typical error bars are indicated in the lower-left. 
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Fig. 11. — Fe II residual vs. Lx/i^FiR- The Fe II residual is defined in the text (see section 4.3). 
The arrows with a circle indicate upper limits for additional lo-BAL QSOs while those without 
a circle are for potential lo-BAL QSOs. For the 14 ROSAT detected objects (filled squares), the 
correlation coefficient is 0.44 with V = 0.117. Typical error bars are indicated in the lower-right. 
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Fig. 12. — IRAS color-color diagram. Our sample galaxies fall in the region between the power- 
law and blackbody lines. The objects close to the blackbody line all have significant [O III] A5007 
blueshifts. 
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Table 1: Sample of IR QSOs Selected from ULIRGs. 
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Note. — The Prefix of the object name indicates the origin of IRAS fluxes. 'F' refers to the IRAS Faint Source 

Catalogue and 'Z' moans the Faint Source Reject File (sec Moshir ct al. 1992). For IR06269-0543, the IRAS fluxes 
come from the IRAS Point Source Catalogue. Units of right ascension are hours, minutes, and seconds, and units 
of declination are degrees, arcminutes, and arcseconds (J2000.0). Col. (4) & (5): far- infrared & infrared luminosity, 
calculated following Sanders & Mirabel (1996); Col. (6): Soft X-ray luminosity (0.2-2.4 keV) normalized to far-infrared 
luminosity; Col.(7); redshift, taken from the references. 



References. — (1) Lawrence et al. 1999; (2) Kim & Sanders 1998; (3) Moran et al. 1996; (4) Zheng et al. 1999. 
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Table 2: Journal of Observations. 







Exp. Time 


Slit 


Seeing 


IRAS Name 


Date 


(second) 


(arcsec) 


(arcsec) 


F00275-2859 


1999 Nov 07 


3600 


3.0 


3.0 


F01572+0009 


1996 Nov 17 


3600 


3.0 


3.0 


F02054+0835 


1998 Oct 22 


2700 


3.0 


2.1 


F02065+4705 


1999 Feb 20 


1800 


2.2 


1.5 


F04415+1215 


1998 Oct 18 


3600 


3.0 


3.5 


IR06269-0543 


1998 Oct 23 


1800 


3.0 


2.1 


F07599+6508 


1997 Mar 16 


3600 


3.0 


1.5 


F09427+1929 


1998 Dec 20 


3600 


2.5 


1.5 


F10026+4347 


1999 Feb 22 


3600 


2.2 


1.5 


F11119+3257 


1999 Feb 21 


2400 


2.2 


1.5 


Z11598-0112 


1997 Mar 12 


5400 


3.0 


3.0 




iqqq Fob 99 


9700 

Li I \J\J 


9 9 

Li. Li 


1 5 


F12265+0219 


1995 Mar 13 


600 


3.0 


3.0 


F12540+5708 


1997 Mar 16 


1200 


3.0 


1.5 


F13218+0552 


1999 Feb 21 


3600 


2.2 


2.0 


F13342+3932 


1999 Feb 22 


1500 


2.2 


1.5 


F15069+1808 


1999 Feb 22 


1200 


2.2 


1.5 


F15462-0450 


1997 Apr 11 


3600 


3.0 


2.0 


F16136+6550 


1998 Oct 20 


1800 


3.0 


2.1 


F18216+6419 


1998 Oct 17 


2700 


2.5 


3.5 


F20036-1547 


1997 Oct 03 


2400 


2.5 


3.5 


F20520-2329 


1998 Oct 22 


2700 


3.0 


2.1 


F21219-1757 


1998 Oct 20 


2700 


3.0 


2.1 


F22454-1744 


1998 Oct 20 


2700 


3.0 


2.1 


F23411+0228 


1997 Oct 04 


4800 


2.5 


3.5 
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Table 3: Profile properties of H/3 and [O III] A5007. 











[0 III]A5007 




FWHM 

J. VV ±±1V± 


"Rlnpcln^fj- 




FWHM 


Blueshift^ 


TRAS Name 




(\rrYl G ■'■^ 




(kms"''^) 


(kms"-*^) 








(A) 


(5) 


(6) 


F00275-2859 


1640 


750 


0.057 


710 




F01572+0009 


2100 


-50 


-0.028 


680 


510 


F02054+0835 


2500; 


1410; 


0.104 


1300; 




F02065+4705 


1410 


450: 


-0.033 


500 




F04415+1215 


1810 


510 


0.180 


1240 




IR06269-0543 


1550 


80 


0.043 


1210 


550 


F07599+6508 


3150 


2030 


0.254 






F09427+1929 


2100 


1640 


0.234 






F10026+4347 


2390 


470 


0.095 


810 




F11119+3257 


1980: 


260: 


0.066 


1480: 


950; 


Z11598-0112 


820 


970 


0.260 


550 




F12134+5459 


740 


280 


0.177 


780 




F12265+0219^ 


3260 


-180: 


-0.004 


1330: 


1020 


F12540+5708 


3130 


1150 


0.101 






F13218+0552 


1180^= 






1530; 


500; 


F13342+3932 


980 


920 


0.446 


540 




r ioUby+ioUo 


looO 


1 on 
— izU 


— O.OoO 


630 




F15462-0450 


1460 


290: 


0.273 


1720: 


1110: 


F16136+6550b 


7270 


-730: 


0.017 


920 


780 


F18216+6419'' 


5190 


-3100; 


-0.293 


1360 




F20036-1547 


1570: 


400: 


0.116 






F20520-2329 


1660 


660 


0.263 


300 


560: 


F21219-1757 


2080 


-40 


-0.045 


1360 


460: 


F22454-1744 


880 


30 


0.028 


790 




F23 11 1+0228 


9 TO" 






1390: 


500; 



"For sources that have an obvious asymmetric [O III] AA4959,5007 profile. 

''Cases with broad H/3 Une; they could be fitted satisfactorily with three Gaussian components. 

'^FWHM of Ha. 

Note. — The uncertainty on the FWHM is typically of order 10% (~200 kms~^). Colons (;) indicate values 
with a relative uncertainty of 20%. Semicolons (;) indicate values with a relative uncertainty 30%-40%. The typical 
uncertainty of the blueshift is < 150 kms^^. Colons (:) indicate values with a relative uncertainty of 30%. Semicolons 
(;) indicate values with a relative uncertainty of larger than 50%. 
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Table 4: Emission Lines Properties. 









Fell 


A5007 


[OIII] 
H/3 


[Oil] 
H/3 


[Nil] 
Ha 


Ha 


H/3 


FeH 


[OHI] 


[OH] 


IRAS Name 


EfB— V) 


Fell 
H/3 


residucil 


PeaJj 


EW 


EW 

(A) 


EW 


EW 


EW 




(■^) 






\^) 




\' ) 


\°) 










(16) 


F00275-2859 


0.41 


1.47 


-0.20 


0.38 


0.14 


0.23 


<0.01 


453.0 


75.6 


110.8 


10.7 


7.5 


F01572+0009 


0.00 


1.13 


0.0 


2.57 


1.05 


0.11 


<0.01 




43.5 


49.0 


45.8 


2.8 


F02054+0835 


0.58: 


2.42: 


-0.06 


0.63: 


0.28: 


<0.16 


0.04: 


169.3: 


33.4: 


80.9: 


9.2: 


<3.1 


F02065+4705 


0.25 


1.65 


0.12 


0.83 


0.35 


<0.06 


0.01 


189.7 


38.4 


63.5 


13.6 


<1.4 


F04415+1215 


0.11 


1.74 


0.0 


0.72 


0.55 


<0.05 


0.12 


209.7 


50.9 


88.3 


28.0 


<1.6 


IR06269-0542 


0.39 


0.57 


0.0 


1.54 


1.02 


0.21 


<0.01 


404.0 


83.6 


47.6 


85.2 


12.2 


F07599+6508 


0.00 


2.75 


-0.21 


<0.08 


<0.02 


<0.04 


<0.01 


258.3 


47.2 


129.6 


<1.1 


<1.3 


F09427+1929 


0.00 


2.98 


-0.09 


<0.08 


<0.02 


<0.02 


0.08 


222.0 


52.0 


155.1 


< 1.1 


<0.6 


F10026+4347 


0.00 


2.07 


0.14 


0.19 


0.10 


<0.01 


<0.01 


149.4 


36.2 


75.0 


3.8 


<0.1 


F11119+3257 


1.08 


1.12 


0.0 


1.01 


0.75 


0.64 


0.12 


181.8 


51.1 


57.2 


,38.5 


25.9 


Z11598-0112 


0.01 


1.71 


0.81 


1.10 


0.75 


0.88 


0.26 


69.8 


13.9 


24.0 


10.5 


8.6 


F12134+5459 


0.14 


1.28 


0.0 


1.13 


0.68 


0.12 


0.17 


148.4 


33.4 


42.7 


22.8 


7.6 


F12265+0219 


0.00 


0.66 


0.0 


0.21 


0.06 


<0.02 


0.07 


299.8 


78.3 


51.4 


4.9 


<0.8 


F12540+5708 


0.65 


1.83 


-0.25 


<0.07 


<0.01 


<0.13 


0.04 


257.2 


35.0 


60.0 


<4.1 


<3.2 


F13218+0552 


0.54 


0.48'' 




1.62: 


2.40=- 




0.28: 


162.7 






47.7: 




F13342+3932 


0.76 


0.73 


0.0 


3.67 


1.09 


0.83 


0.11 


203.2 


35.0 


25.6 


38.0 


17.5 


F15069+1808 


0.01 


0.92 


0.0 


2.04 


1.03 


0.44 


0.18 


189.0 


46.0 


42.2 


47.2 


13.7 


F 15462 -0450 


0.28 


1.32 


0.0 


0.73 


0.66 


0.48 


0.30 


126.7 


48.7 


64.3 


32.2 


13.5 


F16136+6550 


0.17 


0.43 


0.0 


1.52 


0.20 


<0.03 


<0.01 


325.3 


60.0 


26.0 


11.7 


<1.1 


F18216+6419 


0.00 


0.39 


0.0 


1.63 


0.29 


0.04 


0.04 


435.6 


79.9 


30.9 


22.9 


1.7 


F20036-1547 


0.20 


2.74 


-0.09 


<0.10 


<0.01 


<0.15 


0.19 


122.9 


24.2 


66.4 


<0.3 


<2.4 


F20520-2329 


0.46 


2.02 


0.10 


0.72 


0.37 


0.12 


0.33 


103.2 


28.0 


56.5 


10.3 


2.3 


F21219-1757 


0.23 


1.82 


0.05 


0.49 


0.25 


<0.07 


0.32 


153.8 


44.2 


80.5 


11.2 


<1.7 


F22454-1744 


0.21 


1.03 


0.36 


1.52 


1.08 


0.71 


0.17 


166.4 


33.2 


34.3 


35.8 


15.3 


F23411+0228 


0.71: 












0.08 


59.0 











"from Remillard et al. 1993. 

Note. — The typical uncertainty of the emission hne fluxes is about 10-20%. The values flagged with colon 
have uncertainties larger than 40%. Upper limits are given if the emission lines are not convincingly detected. 
No measurement is listed for the emission lines contaminated by noise. Col.(l): source name, same as in table 1; 
Col. (2): color excess; Col. (3): strength of Fe IIA4570 relative to H/9; Col. (4): relative strength of Fe II residual, 
positive indicates excess of Fe 1137,38 and negative indicates excess of Fe 1148,49 (cf Fig. 2); Col. (5): peak height 
of [O III]A5007 relative to that of H/3; Col.(6): strength of [O III] A5007 relative to H/3; Col.(7): strength of 
[O II] A3727 relative to H/3; Col.(8); strength of [N II]A6583 relative to Hq; Co1.(9): equivalent width of Hq; Col.(lO): 
equivalent width of H/3; Col.(ll): equivalent width of Fe IIA4570; Col.(12): equivalent width of [O III] A5007; 
Col.(13): equivalent width of [O II] A3727. 



